k ■ p subband structure of the LaAlOa/SrTiOa interface 
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Heterostructures made of transition metal oxides are new tailor-made materials which are at- 
tracting much attention. We have constructed a 6-band k ■ p Hamiltonian and used it within the 
envelope function method to calculate the subband structure of the LaAlOs/SrTiOa interface, mod- 
eled by an electric field that we treat as a parameter. By use of density functional calculations, we 
determine the k ■ p parameters describing the conduction band edge of SrTiOs: the three effective 
mass parameters, L = 0.6104 eVA 2 , M = 9.73 eVA 2 , TV = -1.616 eVA 2 , the spin orbit splitting 
Aso — 28.5 meV and the low temperature tetragonal distortion energy splitting At = 2.1 meV. We 
find strongly anisotropic non-parabolic subbands. Through quasiclassical quantization for a given 
density, this model allows a direct comparison to the frequency of the Shubnikov-de Haas quantum 
oscillations recently observed in high magnetic fields. Our results are in very good agreement with 
the experiments for an electric field strength F = 0.1 meV/A and a charge density of 7.2 x 10 12 cm~ 2 . 

PACS numbers: 73.20.-r. 71.15.-m,71.20.-b,75.47.-m 
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I. INTRODUCTION 

Oxides of the transition metals, like LaA103 and 
SrTiOs, are intriguing and useful materials, displaying 
many properties attributed to electron correlations ef- 
fects like colossal magnetoresistance, ferro- and antifcrro- 
magnetism, metal-insulator transitions, and high and low 
temperature superconductivity. Recently, pulsed laser 
deposition growth methods succeeded to combine differ- 
ent oxides in heterostructures with layers as thin as 10- 
20 lattice parameters and relatively sharp interfaces. 1 ' 2 
These samples resemble the more familiar semiconduc- 
tor heterostructures since the different band gaps of the 
two materials and their band line-up at the interface 3 
can lead to quantization of the electronic states into two- 
dimensional levels, opening the way to the typical two 
dimensional electron systems phenomena like the quan- 
tum Hall effect and conductance quantization. 

Much interest in these structures has been triggered by 
the observation in magnetotransport of a high-mobility 
electron gas at the LaAlOs/SrTiOs interface, 2,4 an un- 
expected feature in view of the fact that the constituent 
materials are both insulators. Several models have been 
proposed to justify the presence of this charge distribu- 
tion at the interface, leading to widely spread estimates of 
the resulting subband energy spacing (from few to hun- 
dreds of meV) and either atomic or extended states at 
the interface. 3-8 As yet, the origin of the charge carriers 
is still unclear. 

In this paper we present band structure calculations 
with a 6-band k ■ p approach for the conduction band 
of the bulk materials as proposed in Ref. 9 and match- 
ing of the envelope functions at the interface for the 
LaAlOs/SrTiOs heterostructures, an approach success- 
fully used for semiconductors. In this effective mass 
approach, which is very convenient in view of the very 
large number of atoms in the unit cell, one writes an 
effective Schrodingcr equation for the envelope of the 
atomic states in the heterostructure. The SrTiOs k ■ p 



parameters are determined by fitting the bulk band struc- 
tures calculated within density functional theory, yielding 
three effective mass parameters L, M, N, the spin orbit 
splitting Aso and the low temperature tetragonal dis- 
tortion energy splitting Ay. The k ■ p envelope function 
method not only allows us to determine the quantized 
energy levels but also to calculate the in-plane dispersion 
determining the density of states and the quasiclassical 
quantization in perpendicular magnetic fields. 

Our calculations are particularly relevant for the analy- 
sis of recent magnetotransport data on the high-mobility 
electron gas at the LaAlOs/SrTiOs interface 10 which 
shows multisubband Shubnikov-de Haas (SdH) oscilla- 
tions with a few meV subband separation. We find 
anisotropic non-parabolic bands with quite different in- 
plane masses and an energy dependent density of states 
that is important for the interpretation of the SdH ex- 
periments, which usually assumes parabolic and isotropic 
bands. Moreover, as the relevant energy spacings at the 
interface are meV, incorporation of spin-orbit (SO) cou- 
pling is crucially important. We assume a positive charge 
localized at the interface giving rise to an electric field in 
the SrTiOs. The strength of this electric field is taken 
as a constant parameter since we do not incorporate any 
lattice and electron screening in our model. For charge 
density n = 7.2 x 10 12 cm -2 and electric field strength 
F = 0.1 meV/A we can reproduce the magnetotransport 
data quite precisely using only 5 relevant k p parameters. 

In section II we model the bulk conduction band struc- 
ture around the T point using the fe • p approach. In 
section III we find the k ■ p parameters by fitting the 
bulk bands calculated within density functional theory 
(DFT). We discuss in detail the importance of spin or- 
bit coupling. In section IV we use the envelope func- 
tion method to calculate the subband energy structure 
of SrTiOs quantum wells. In section V we introduce an 
electric field to account for the polar structure of the in- 
terface and compare our results to the SdH experiments. 
Section VI presents a summary and discussion of this 
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work. 



II. BULK BAND STRUCTURE MODEL 

We construct a model of the band structure of the 
LaAlOs/SrTiOa interface around the T point based on 
the k p approach. This requires a description of the band 
edges of the constituents and their offset. By appropriate 
matching of the wave functions at the interface, one can 
calculate the energy levels. Both SrTiOs and LaAlC>3 are 
insulators, with an experimental band gap of 3.2 and 5.6 
eV respectively. 11 ' 12 In Fig. 1 we show the energy bands 
of the two constituents, calculated as described in sec- 
tion III, with the symmetry and atomic character at T. 
Based on the position and symmetry of the bands in the 
two materials, we argue that a description of the SrTiC>3 
band edge only is sufficient. In fact, the conduction band 
of SrTi03 has a minimum at the T point constituted by 
the 6- fold degenerate (including spin) Ti ti g bands. The 
tig bands are a subset of the d bands, d xy ,d yz and d zx , 
split by the crystal field from the other d bands, the e g 
states which lie ~ 2 eV higher. In the cubic perovskite 
structure the ti g bands at T transform as rj 5 . 13 These 
bands should be matched to states of the same symme- 
try in the LaAlC>3 layer. The lowest in the LaAlC>3 
are located well above the band gap and since the band- 
offset is expected to be of type I, 3 even further above the 
tig states in SrTiC>3. Moreover, in LaAlC>3 the states of 
this symmetry are localized on La, namely on a different 
crystalline position than the Ti in the SrTi03 , with con- 
sequent small overlap. Therefore, we can safely neglect a 
spreading of the Ti ti g states in the LaA103 and assume 
an infinite potential barrier at the interface. This greatly 
simplifies the problem, avoiding explicit matching of the 
wave functions in the two materials. 

From now on we focus on the 6-fold degenerate ti g 
states in SrTiC>3 which have symmetry like the p- 
states in III-V semiconductors with zincblende struc- 
ture 13 . Note that, due to the crystal field, the orbital 
momentum I = 2 of the d-states is no longer a good 
quantum number. The spin-orbit splitting will be con- 
sidered as a perturbation on this level structure and will 
split the 6- fold degeneracy into a 4- fold J = 3/2 mul- 
tiplet of symmetry and a 2-fold J = 1/2 multiplct 
of symmetry , exactly as in the valence band of III-V 
semiconductors . 13 

Following Ref. 9 we describe the bulk Ti ti g bands 
around T by means of a k ■ p Hamiltonian depending on 
3 effective mass parameters L, M, N, the spin-orbit split- 
ting Aso an d the low temperature tetragonal distortion 



energy splitting 

(L,M,N) + H so (Aso) + H T (A T ). (1) 

We choose as basis functions the six k = ti g states 
\X t),\Y t),\Z t),\X i),\Y i),\Zi) corresponding to 
d yz ,d zx and d xy with both spin up and spin down, al- 
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Figure 1. Bulk band structure of (top) cubic SrTiC>3 and (bot- 
tom) cubic LaA103 calculated without spin-orbit interaction. 
The states at T are labeled by the symmetry representation 
and by the dominant atomic character. The La f states have 
symmetry , r^ 5 , r^" 5 with increasing energy. Notice that 
only the F^ 5 Ti tig bands are considered to calculate the en- 
ergy subbands at the LaAlOa/SrTiOa interface (see text). 



though H cu bi C and Ht do not depend on spin. 
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where H T + describes the r^~ 5 bands in a cubic crystal 
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In this Hamiltonian the interaction with remote bands 
is taken into account by the effective mass parameters 
L, M and TV. TV quantifies the coupling between the ti g 
states via states with other symmetries. 

At T = 110 K SrTiC>3 undergoes a structural phase 
transition to a tetragonal symmetry. At T — 4.2 K the 
TiO octahedra have rotated 2.1° around the tetragonal 
axis. 14 We include this distortion in our model because 
we want to compare our results to low T transport exper- 
iments. For the heterostructures, we choose the tetrag- 
onal axis along the growth axis which we call the z axis 
throughout this article. Ht has only two non-zero ma- 
trix elements (Z t \V\Z f) = {Z 1 \V\Z \) = A T that, 
at r, shift the \Z) band A T above \X) and \Y). 

The effect of spin-orbit interaction on the t2 g bands 
is described only by the splitting A S o Ht T between the 
4-fold T+ and 2-fold T+ level: 



explicitly the matrix H in the \ Jrrij) basis of Eq.(5): 
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which lowers the J = 3/2 multiplet by Ago/3 and raises 
the J = 1/2 multiplet by 2Ago/3. In the following, for 
simplicity, we take the origin of energy at the J = 3/2 
multiplet which is the minimum of the conduction band. 
Since the lowest states are almost purely eigenstates of 
the total angular momentum J, it is convenient to write 
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r =-(L + 2M){k 2 x + k 2 + ki) + A so + -A T . 

The effective mass parameters L, M and N and the en- 
ergy splittings At and A$o have to be determined by 
experiments or calculations. As discussed in detail in 
Ref. 9 the experimental data for SrTiC>3 is rather scarce 
and contradictory and therefore we have performed DFT 
calculations which also allow to study each term in the 
Hamiltonian separately. 



III. BAND STRUCTURE PARAMETERS 
A. DFT calculations 

We perform first-principles calculations in the frame- 
work of DFT (Refs. 15 and 16) employing the Perdew- 
Burkc-Ernzcrhof (PBE) generalized gradient approxima- 
tion (GGA). 17 The projector augmented wave (PAW) 
method 18 ' 19 is applied, as implemented in the Vienna 
ab initio simulation program (VASP). 20 ' 21 We use stan- 
dard PAW data sets as provided with the VASP package, 
which have for Sr a frozen [Ar]3d 10 core, for Ti and Al a 
frozen [Ne] core, and for La a frozen [Kr]4d 10 core. The 
La data set includes /-channel augmentation. Calcula- 
tions with a harder data set (smaller PAW core radii) 
for oxygen, confirm our results. VASP uses spinors for 
calculations with spin-orbit coupling in the Kohn-Sham 
(KS) Hamiltonian. 

A kinetic energy cutoff of 400 eV is employed for the 
plane wave expansion of the KS orbitals. The calcula- 
tions are done with a T-centered 12 x 12 x 12 k-mesh. 22 
We use experimental lattice parameters at room temper- 
ature: a = 3.905 A for SrTi0 3 (Ref. 23) and a = 3.789 A 
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k ■ p parameters 


numerical value 


calculation 




n fil ftAC?) o\lk z 
\j . u i-U^±y^ j cv -r\- 


piiKi/i via S[ 1 


M 


9.73(2) eVA 2 


cubic no SO 


N 


-1.616(4) eVA 2 


cubic no SO 


Aso 


28.5 meV 


cubic with SO 


At 


2.1 meV 


tetragonal no SO 



Table I. kp parameters found by fitting the DFT band struc- 
ture of SrTiOa. 



for LaA103. 24 For body-centered tetragonal SrTiOs we 
use a = 3.896 A and c = 3.899 A. 23 

To find the character of the bands, we project the KS 
states onto spherical harmonics, within a radius r p . For 
SrTiC>3 we choose r p w 1.2 A for all elements. 

The resulting band structures for cubic SrTiOs and 
cubic LaA103 are shown in Fig. 1. 

B. Fit effective mass parameters 

The energy band dispersion given by Eq. 3 along the 
three high symmetry lines k = (k,0,0), fen = (k, k, 0) 
and k = (k, k, k) can be written analytically in terms of 
L, M and N. These effective mass parameters are then 
determined by fitting to the DFT band structure of cubic 
SrTiOs without spin-orbit coupling. The k x -, k y - and /in- 
directions are equivalent because of the cubic symmetry. 
The resulting values for L, M and TV are listed in Table I. 
In Fig. 2a the DFT bands arc compared to the kp model 
along k z and k\ \ . Along k z there is one very flat band with 
to* = 6.24 mo- This unusual heavy mass originates from 
the fact that the \Z) state lies in the xy-plane and does 
not extend along the z direction, as shown in the inset 
of Fig. 2a. Along fc|| we find three distinct dispersions 
given by 1/2(L + M + N)k 2 , \/2(L + M-N)k 2 and Mk 2 
which indicate that N ^ 0, in contrast to what has been 
used in Ref. 8. 



C. Calculation of Aso and At 

In this Section we discuss how the inclusion of SO 
coupling and of a tetragonal distortion modify the band 
structure of Fig. 2a. To determine Aso we calculate 
the DFT energy band structure of cubic SrTiC>3 includ- 
ing spin-orbit coupling. We find Aso = 28.5 meV, in 
agreement with other calculations. 25 Fig. 2b shows the 
resulting 4-fold J = 3/2 states and the 2- fold J = 1/2 
split-off bands at T. Although the latter is located Aso 
above J = 3/2, it affects the dispersion of the |3/2± 1/2) 
band to which it is coupled at finite k z . Along k z , the 
|3/2 ± 3/2) and the |3/2 ± 1/2) arc the light and heavy 
mass electrons respectively. This is opposite to the well- 
studied case of p-type valence bands in III-V semicon- 
ductors, which can be easily understood by noting that 
while the p z lobe extends mainly in the z direction, the 



\Z) state extends only in the zy-planc. Notice that the 
spin-orbit interaction makes the effective heavy electron 
mass much lighter, bringing it from m* = 6.24 mo to 
to* = 1.17 mo and with an almost linear dispersion at 
larger k z . 

The tetragonal distortion energy splitting At is deter- 
mined by calculating the band structure of the tetrag- 
onally distorted SrTiOs without spin-orbit interaction. 
We find At = 2.1 meV, close to the values found in 
Refs. 8 and 26. The tetragonal distortion breaks the 
symmetry between k x and k z as is shown for the case 
of tetragonally distorted SrTiOs including SO in Fig. 2c. 

In summary, the k ■ p model is in excellent agreement 
with the DFT calculations at least up to k ~ 0.1 A -1 , 
which is the relevant range for confinement over lengths 
larger than ~ 40 A. 



IV. QUANTUM WELL SUBBANDS 

In the previous section we have shown that the k ■ p 
method describes very well the conduction band of bulk 
SrTiOs. A heterostructure leads to a potential profile 
along the growth direction V(z), leading to quantization 
of k z . The wave function can then be written as 

^ fc|1)Jk ,(r)=e 4k ii-'iitt(r).^) (9) 

with k|| = (k x ,k y ,0), jpi I = (x,y,0) and the basis func- 
tions Ui(r) have the periodicity of the bulk unit cell. 
The components of 4>{z) 1 <\>\ (z) . . . 4>6(z) are the enve- 
lope functions replacing e lkzZ which vary slowly at the 
scale of the unit cell. They can be found by solving the 
eigenvalue equation: 

\h (k z -> -i^j + V(z)l} 4>(z) = E(k)4>(z) (10) 

From now on we will always use the full Hamiltonian 
H, Eq. (1). We solve this equation by the finite difference 
method, where we discretize the envelope wave function 
on an equispaced grid in real space. 27 

For illustration, in Fig. 3 we show the quantized energy 
levels in a 100 A wide SrTiOs quantum well (QW) with 
infinite barriers. With increasing energy, the first, third 
and fourth subbands derive from the heavy |3/2 ± 1/2) 
bands and have a mixed |3/2±l/2) and |l/2±l/2) char- 
acter, whereas the second subband has pure 1 3/2 ± 3/2) 
character. In the in-plane direction, the heavy electrons 
have a light effective mass and vice versa. This leads 
to strongly non-parabolic dispersion of the in-plane sub- 
bands. In particular we see a strong avoided crossing 
between the first two subbands at finite in-plane wave 
vectors. As a consequence, the equal energy contours are 
strongly anisotropic and energy dependent, a property 
which is important for the interpretation of SdH experi- 
ments that we will discuss in the next section. 
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Figure 2. Dispersions of the t2 g bands in cubic SrTi03 according to DFT calculations (dots) and to the k ■ p model (lines) 
with the parameters given in Table I. (a) cubic SrTi03 without SO (H = H cu bi c ) along k z and fcii = (k,k,0). Along k z the 
steepest band is 4- fold degenerate. The inset shows the d yz orbitals (upper) and d xy orbitals (lower). The three distinct bands 
along fc|| indicate N / 0. (b) cubic SrTiC>3 with SO (H — H cu bi C + Hso)- Note that the bands are shifted by Aso/3 so that 
the j = 3/2 multiplet edge is at E — 0. (c) tetragonally distorted SrTiOs with SO (H = H cu bi c + Hso + H& T ). Note that k x 
and k z are no longer equivalent. 
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Figure 3. Left: quantized levels (horizontal lines) and corre- 
sponding k = envelope functions of a 100 A SrTiOa QW. 
The second subband (dashed lines) has a pure |3/2 ± 3/2) 
character. The envelope function of the states derived from 
the |3/2 ± 1/2) band (solid lines) are mixed and we plot 
both 1 3/2 ± 1/2) and 1 1/2 ± 1/2) components: the modulus 
squared of the |3/2 ± 1/2) component is 98%, 88% and 80% 
for increasing energy. Right: in-plane dispersion. Notice the 
non-parabolic dispersion due to the strong avoided crossing 
between the first and second subbands.The horizontal line at 
E = 16 meV corresponds to the energy for which the, strongly 
anisotropic, equal energy contours are shown in the inset. 



V. EFFECT OF ELECTRIC FIELD AND 
COMPARISON TO MAGNETOTRANSPORT 
EXPERIMENTS 

The LaAlOa/SrTiOs interface leads to a potential pro- 
file V(z). We assume a uniform positive charge density 
at the interface which gives rise to a triangular potential 



in the SrTiC>3. As already argued, the interface can be 
modeled by an infinite barrier. 

V(z) =Fz if z>0 

V(z) = 00 if z < ^ ' 

The strength of the electric field is related to the charge 
density at the interface which is not directly experimen- 
tally accessible. Moreover, screening of the field by ions 
and electrons is complex and lies beyond the scope of our 
model. We therefore treat the electric field strength F as 
a constant parameter in our calculations. 
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Figure 4. Dependence of the first 10 quantized levels on F. 
Dashed lines: |3/2 ± 3/2) states. The other subbands have 
a mixed |3/2 ± 1/2) and 1/2 ± 1/2) character. Notice the 
meV spacing in this range of fields. The vertical line is at the 
field F = 0.1 meV/A for which the subband dispersion is 
presented in Fig. 5. 

In Fig. 4 we show the dependence of the first 10 quan- 
tized levels on the electric field strength F which repre- 
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scnts the slope of the confining triangular well. We report 
with dashed lines the pure |3/2± 3/2) states. The bands 
originating from the 1 3/2 ± 1/2) bulk bands (solid lines) 
have a mixed |3/2± 1/2) and |l/2± 1/2) character. The 
modulus squared of the 1 1/2 ± 1/2) is below 10% for all 
subbands in the figure. Notice that, in this range of elec- 
tric fields, the spacing of the levels is of the order of meV 
as found by tight binding in Rcf. 8 in the low density 
regime with n of the order of 10 12 — 10 13 cm -2 . We see 
that the dependence is rather complex with level cross- 
ings due to the spin-orbit interaction and the tetragonal 
distortion. In Fig. 5 we show the envelope functions and 
subband dispersion calculated for electric field strength 
F = 0.1 meV/A, because, as we will discuss later, this 
value gives the best agreement with the SdH experiments. 
We see that the envelope functions at k = have the 
asymmetric shape corresponding to the Airy functions. 
For this value of the field, the extension of the first sub- 
bands is of the order of 100 A which justifies the choice of 
the envelope function method. The in-plane dispersion, 
as for a QW, is strongly non-parabolic and anisotropic, 
particularly due to the strong avoided crossing between 
the first and second subband at finite k. Notice that, 
as already found in the bulk, spin-orbit coupling yields 
an almost linear dispersion of the lowest subband which 
gives rise to a heavy effective mass. 

An additional feature is represented by the small split- 
tings which result from spin-orbit coupling in a potential 
without inversion symmetry that breaks the Kramer's 
degeneracy. 28 Notice that the splittings disappear for 
N = because in that case the Hamiltonian of Eq. (5) 
splits into two equivalent 3x3 blocks. The spin splitting 
has not been found in Ref. 8 where the mixing parameter 
N is neglected. 

SdH oscillations are periodic in the inverse magnetic 
field. Using quasiclassical quantization one can relate 
the frequencies / measured in an experiment directly to 
the Fermi surface S in fc-space 29 



27re/ 



(12) 



We use this relation to compare our calculations to 
the experimental results. 10 We calculate the surfaces and 
corresponding frequencies for various values of F as a 
function of energy. We average over the spin-split bands, 
as these small splittings cannot be resolved in experi- 
ments. In Fig. 6 we show the calculated frequencies as 
a function of the chosen Fermi energy for F = 0.1 meV, 
together with the frequencies that were measured for one 
of the three samples in Ref. 10 (sample 2). We see that 
at this value of the electric field, at Ep = 18.6 meV our 
calculations are in very good agreement with the experi- 
mental values. The electron density for this Fermi energy 
is n = 7.2 x 10 12 cm" 2 . 

From the temperature dependence of the SdH oscilla- 
tions one can extract an average effective mass at the 
Fermi energy. As the bands are neither parabolic nor 
isotropic, the average effective mass at the Fermi energy 
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Figure 5. Left: quantized levels (horizontal lines) and cor- 
responding k — envelope functions of a triangular SrTi03 
well with slope F — 0.1 meV/A indicated by a solid line. 
The second and fifth subbands (dashed lines) have a pure 
3/2 ±3/2) character. The envelope function of the states de- 
rived from the 1 3/2 ± 1/2) band (solid lines) are mixed and we 
plot both 1 3/2 ± 1/2) and 1/2 ± 1/2) components: the mod- 
ulus squared of the |3/2 ± 1/2) component is 98%, 96% and 
94% for increasing energy. Right: in-plane dispersion. Notice 
the almost linear dispersion of the first subband away from 
r. A small spin splitting is due to the non symmetric confin- 
ing potential (see text). The horizontal line at E — 18.6 meV 
corresponds to the energy for which the equal energy contours 
are shown in the inset. 



is difficult to compute for the subbands we have calcu- 
lated. One way to proceed is to calculate the density of 
states (DOS) and relate this to the effective mass. For 
a two-dimensional system with parabolic isotropic bands 
the DOS is constant 



DOS2D = 



(13) 



explicitly counting each spin. We calculate the DOS for 
each subband by use of the k space energies on a fine 
grid, and building a normalized Gaussian with 0.8 meV 
width around each point. The resulting DOS is not con- 
stant, as the subbands are neither parabolic nor isotropic. 
The DOS represents an average of the subband dispersion 
over all k points with the same energy, and it can there- 
fore be related to the energy dependent effective mass by 
inverting Eq. (13). In this way, we can determine for each 
subband the average effective mass at the Fermi energy, 
and compare it to the values that were found in Rcf. 10. 

In Fig. 7 we show the DOS and the corresponding ef- 
fective mass as a function of energy. The precise values of 
the effective masses are very sensitive to the chosen Ep 
but the order of heavy and light masses is a robust fea- 
ture. In Ref. 10 the authors find masses corresponding to 
the four frequencies reported in Fig. 6 of 2.0 mo, 0.9 mo, 
0.9 mo, 0.9 mo with decreasing frequency. We also find 
that the mass of the lowest subband is more than twice 
that of the following three subbands. For Ef = 18.6 meV 
and averaging over the spin split states we find 1.5 mo, 
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0.5 mo, 0.5 mo, 0.5 mo which gives a satisfactory agree- 
ment. Note that the spin-split bands can have quite dif- 
ferent effective masses, for example the effective masses 
of the third spin split subband are m* = 0.42 mo and 
m* = 0.52 mo. 



100 




12 14 16 15 
E F (meV) 

Figure 6. Dependence of the Fermi surface (right axis) or of 
the frequency of quantum oscillations (left axis) on the Fermi 
energy for all occupied subbands calculated with F = 0.1 
meV/A. The Fermi surface has been calculated by averaging 
over the small spin splittings (see Fig. 5) . The horizontal 
lines correspond to the frequencies observed in SdH exper- 
iments of Ref. 10 for sample 2: fi — 83 Tesla, fa = 36 
Tesla, / 3 = 20 Tesla, f 4 = 8 Tesla. Notice that at the en- 
ergy Ef = 18.6 meV, indicated by a vertical solid line, the 
first four frequency coincide with the measured values. 



VI. SUMMARY AND DISCUSSION 



possibly due to defects, exists at the interface, creating 
an electric field that confines conduction electrons in the 
SrTiOa. The parameters necessary to describe the bulk 
have been determined by fitting to the band structure 
calculated by DFT. We find several occupied subbands 
a few meV apart with different masses as also found ex- 
perimentally. We even find an excellent agreement with 
one of the sample measured. For these new oxide het- 
erostructures reproducibility between samples is not yet 
very good. In all measured cases, however, the subband 
spacing is of the order of meV, which supports the phys- 




io 12 14 16 18 20 

E F (meV) 

Figure 7. Dependence of the DOS (right axis) and related 
average effective masses (left axis) on the Fermi energy for 
F — 0.1 meV/A. For reference, the dashed line gives the 
constant density of states of a free electron (m* = mo) in 
two dimensions. Notice that the DOS deviates from the free- 
electron like behavior, particularly for the lowest spin-split 
subband that grows almost linearly with Ef- The vertical 
dotted line is drawn at Ef = 18.6 meV. 



ical picture of an electronic distribution at the interface 
spread over many atomic layers. The best agreement is 
found for a relatively low value of the electric field. We 
hope that our results will stimulate further research on 
this topic. 



In summary, we have calculated the subband struc- 
ture at the LaAlOa/SrTiOa interface with the k p enve- 
lope function method. We have calculated the frequency 
of SdH oscillations and compared the results to recent 
experimental values. We have assumed, based on pre- 
vious literature, that some positive charge distribution, 
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